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Abstract 

We establish explicit duality transformations for systems of M 
g-state Potts models coupled through their local energy density, gen- 
eralising known results for M = 1,2,3. The M-dimensional space 
of coupling constants contains a selfdual sub-manifold of dimension 
Dm = [M/2]. For the case M = 4, the variation of the effective 
central charge along the selfdual surface is investigated by numerical 
transfer matrix techniques. Evidence is given for the existence of a 
family of critical points, corresponding to conformal field theories with 
an extended Sm symmetry algebra. 

For several decades, the g-state Potts model has been used to model 
ferromagnetic materials |IJ, and an impressive number of results are known 
about it, especially in two dimensions 0, [|, [§]. More recently, its random- 
bond counterpart has attracted considerable attention || , primarily because 
it permits one to study how quenched randomness coupling to the local 
energy density can modify the nature of a phase transition. 

But despite the remarkable successes of conformal invariance applied 
to pure two-dimensional systems, the amount of analytical results on the 
random-bond Potts model is rather scarce. Usually the disorder is dealt with 
by introducing M replicas of the original model, with mutual energy-energy 
interactions, and taking the limit M — ■> 0. The price to be paid is however 
that the resulting system loses many of the properties (such as unitarity) 
that lie at the heart of conventional conformal field theory || |7[]. 

Very recently, an alternative approach was suggested by Dotsenko et al 
M. These authors point out that the perturbative renormalisation group 
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M (effectively an expansion around the Ising model in the small parameter 
e = q — 2) predicts the existence of a non-trivial infrared fixed point at 
interlayer coupling oc — ej (M — 2) + 0(e 2 ), so that the regions M < 2 and 
M > 2 are somehow dual upon changing the sign of the coupling const antQ. 
More interestingly, for M = 3 they identify the exact lattice realisation of 
a critical theory with exponents consistent with those of the perturbative 
treatment, and they conjecture that this generalises to any integer M > 
3. Their proposal is then to study this class of coupled models, which are 
now unitary by definition, and only take the limit M — > once the exact 
expressions for the various critical exponents have been worked out. One 
could hope to attack this task by means of extended conformal field theory, 
thus combining the Z q symmetry of the spin variable by a non-abelian Sm 
symmetry upon permuting the replicas. 

Clearly, a first step in this direction is to identify the lattice models corre- 
sponding to this series of critical theories, parametrised by the integer M > 3. 
For M = 3 this was achieved || by working out the duality relations for M 
coupled Potts models on the square lattice, within the M-dimensional space 
of coupling constants giving rise to Sm symmetric interactions amongst the 
lattice energy operators of the replicas. Studying numerically the variation 
of the effective central charge 0] along the resulting selfdual line, using 
a novel and very powerful transfer matrix technique, the critical point was 
unambiguously identified with one of the endpoints of that line. 

Unfortunately it was hard to see how such duality relations could be 
extended to the case of general M. The calculations in Ref. || relied on a 
particular version [|ll|] of the method of lattice Fourier transforms |12|], already 



employed for M = 2 two decades ago [13|. Though perfectly adapted to the 



case of linear combinations of cosinoidal interactions within a single (vector) 
Potts model ||12|| , this approach led to increasingly complicated algebra when 
several coupled models were considered. Moreover, it seemed impossible to 
recast the end results in a reasonably simple form for larger M. 

In the present publication we wish to assess whether such a scenario of a 
unique critical point with an extended Sm symmetry can indeed be expected 

1 The case M — 2 is special: For q — 2 (the Ashkin- Teller model) the coupling presents 
a marginal perturbation, giving rise to a halflinc of critical points along which the critical 
exponents vary continuously On the other hand, for q > 2 where the perturbation 
is relevant, the model is still integrable, but now presents a mass generation leading to 
non-critical behaviour 101. 
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to persist in the general case of M > 3 symmetrically coupled models. We 
explicitly work out the duality transformations for any M, and show that they 
can be stated in a very simple form [Eq. (|9])] after redefining the coupling 
constants. 

The lattice identification of the M = 3 critical point in Ref. || crucially 
relied on the existence of a one-parameter selfdual manifold, permitting only 
two possible directions of the initial flow away from the decoupling fixed 
point. We find in general a richer structure with an [M/2] -dimensional self- 
dual manifold. Nonetheless, from a numerical study of the case M = 4 we 
end up concluding that the uniqueness of the non-trivial fixed point can be 
expected to persist, since the decoupling fixed point acts as a saddlepoint of 
the effective central charge. 

Consider then a system of M identical planar lattices, stacked on top of 
one another. On each lattice site i, and for each layer /x = 1,2,..., M, we 
define a Potts spin that can be in any of q = 2, 3, . . . distinct states. The 
layers interact by means of the reduced hamiltonian 

w = E««> (i) 

(ij) 

where (ij) denotes the set of lattice edges, and an Sm symmetric nearest- 
neighbour interaction is defined as 

M i m 

«o=-E*m E n^(4 w Vj w) )- ( 2 ) 

m=l Mi 7^27^ •••Mm J=l 

By definition the primed summation runs over the rjf) terms for which the 
indices 1 < \i\ < M with I — 1,2, ... ,m are all different, and 5(x, y) — 1 if 
x = y and zero otherwise. 

For M = 1 the model thus defined reduces to the conventional Potts 
model, whilst for M = 2 it is identical to the Ashkin- Teller like model consid- 
ered in Ref. |[L3|| , where the Potts models of either layer are coupled through 



their local energy density. For M > 2, additional multi-energy interactions 
between several layers have been added, since such interactions are gener- 
ated by the duality transformations, as we shall soon see. However, from the 
point of view of conformal field theory these supplementary interactions are 
irrelevant in the continuum limit. The case M = 3 was discussed in Ref. IH1. 
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By means of a generalised Kasteleyn-Fortuin transformation [II] the local 
Boltzmann weights can be recast as 



M i 

ex P (-^,) = n n 

771=1 fJ-l^^^—Hm 

In analogy with the case of M — 1, the products can now be expanded so 
as to transform the original Potts model into its associated random cluster 
model. To this end we note that Eq. (|3|) can be rewritten in the form 

M I m 

Op(-Wy) = b + £ b m £ I] ^\ (4) 

m=l Atl^M2^— Mm i=l 

defining the coefficients {& m }m=o- The latter can be related to the physical 
coupling constants {K m }^ =1 by evaluating Eqs. (|3|) and (|J) in the situation 
where precisely m out of the M distinct Kronecker 5-functions are non-zero. 
Clearly, in this case Eq. (|^) is equal to e Jm , where 

J m = t(f)K k (5) 

for m > 1, and we set Jo = i^o — 0. On the other hand, we find from Eq. (0j) 
that this must be equated to J2T=o This set of M + l equations can be 

solved for the bk by recursion, considering in turn the cases m — 0,1, ... , M. 
After some algebra, the edge weights bk (for k > 0) are then found as 

h= E(-ir +fc (^y™. w 

The partition function in the spin representation 

z = J2ReM-n lJ ) (7) 

can now be transformed into the random cluster representation as follows. 
First, insert Eq. (B) on the right-hand side of the above equation, and imagine 
expanding the product over the lattice edges (ij). To each term in the 
resulting sum we associate an edge colouring Q of the M-fold replicated 



l+(e^-l)n% 



1=1 



(3) 
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lattice, where an edge (ij) in layer m is considered to be coloured (occupied) 
if the term contains the factor 5{a ( ( n \ c r j m ' ) ), and uncoloured (empty) if it 
does not. [In this language, the couplings correspond to the local energy 
density summed over all possible permutations of precisely k simultaneously 
coloured edges.] 

The summation over the spin variables {a} is now trivially performed, 
yielding a factor of q for each connected component (cluster) in the colouring 
graph. Keeping track of the prefactors multiplying the 5-functions, using 
Eq. (f|), we conclude that 

M 

z = J2U i Cm € m , (8) 

Q m=l 

where C m is the number of clusters in the mth layer, and B m is the number 
of occurencies in Q of a situation where precisely m (0 < m < M) edges 
placed on top of one another have been simultaneously coloured. 

It is worth noticing that the random cluster description of the model has 
the advantage that q only enters as a parameter. By analytic continuation 
one can thus give meaning to a non-integer number of states. The price to 
be paid is that the C m are, a priori, non-local quantities. 

In terms of the edge variables b m the duality transformation of the par- 
tition function is easily worked out. For simplicity we shall assume that 
the couplings constants {K m } are identical between all nearest-neighbour 
pairs of spins, the generalisation to an arbitrary inhomogeneous distribu- 
tion of couplings being trivial. By analogy with the case M = 1, a given 
colouring configuration Q is taken to be dual to a colouring configuration 
Q of the dual lattice obtained by applying the following duality rule: Each 
coloured edge intersects an uncoloured dual edge, and vice versa. In partic- 
ular, the demand that the configuration (/ full with all lattice edges coloured 
be dual to the configuration Q cm _ p t y with no coloured (dual) edge fixes the 
constant entering the duality transformation. Indeed, from Eq. (|8]), we find 
that (?f u ii has weight q M bfj, where E is the total number of lattice edges, 
and {?empty is weighted by q MF b F ', where F is the number of faces, includ- 
ing the exterior one. We thus seek for a duality transformation of the form 
q MF b$ Z({b m }) = q M bfj Z ({b m }) , where for any configuration Q the edge 
weights must transform so as to keep the same relative weight between Q 
and (?f u ii as between Q and Q empty 
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An arbitrary colouring configuration Q entering Eq. ([8]) can be generated 
by applying a finite number of changes to Q{ u \h i n which an edge of weight 6m 
is changed into an edge of weight b m for some m = 0,1, . . . , M — 1. By such a 
change, in general, a number k < M — m of pivotal bonds are removed from 
the colouring graph, thus creating k new clusters, and the weight relative 
to that of £?f u ii will change by q k b m jbu- On the other hand, in the dual 
configuration Q a number M — m — k of clusters will be lost, since each of 
the k new clusters mentioned above will be accompanied by the formation 
of a loop in Q. The weight change relative to (/empty therefore amounts to 
bM-m/(boq M ~ m ~ k ). Comparing these two changes we see that the factors of 
q k cancel nicely, and after a change of variables m — > M — m the duality 
transformation takes the simple form 



q m b 



M—m 



for m — 0,1, 



M, 



(9) 



the relation with m = being trivial. 

Selfdual solutions can be found by imposing b r< 



rise to only 



A/+1 



b m . However, this gives 



independent equations 



>M-m 



q 



M/2-mi 

Or, 



for m 



0,1,..., 



M-l 



(10) 



and the M-dimensional parameter space {b m }^ =1 , or {K m }^_-,, thus has a 



m=l ' 

In particular, the ordinary 
2 PH and 



selfdual sub-manifold of dimension Dm 

Potts model (M = 1) has a unique selfdual point, whilst for M 
M = 3 H one has a line of selfdual solutions. 

Our main result is constituted by Eqs. (H) and (H) relating the physi- 
cal coupling constants {K m } to the edge weights {b m }, in conjunction with 
Eqs. (0) and (10) giving the explicit (self) duality relations in terms of the 
latter. 

Since the interaction energies entering Eq. (§) are invariant under a simul- 
taneous shift of all Potts spins, an alternative way of establishing the duality 
transformations procedes by Fourier transformation of the energy gaps [[11] 



This method was used in Refs. [13] and [§] to work out the cases M = 2 and 
M = 3 respectively. However, as M increases this procedure very quickly be- 
comes quite involved. To better appreciate the ease of the present approach, 
let us briefly pause to see how the parametrisations of the selfdual lines for 
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M = 2,3, expressed in terms of the couplings {K m }, can be reproduced in a 
most expedient manner. 

For M = 2, Eq. flTOD gives b 2 = q, where from Eqs. (|5|) and (|5|) b 2 = 

e 2K 1+ K 2 _ 2e K, + L Thug 



,K 2 



2e Kl + 



(11) 



in accordance with Ref. [13 . Similarly, for M = 3 one has b\ = qb 2 /b% = 
b 2 j^fq with b x = e Kl - 1, b 2 as before, and b 3 = e ^i+3K 2 +K 3 _ ^ e 2K 1+ K 2 + 

3e Kl — 1. This immediately leads to the result given in Ref. p|: 



,k 2 



,K 3 



2 + V9)e /fl -(l + V9) 



(12) 



3(e^-l)(l + 7g) + g 3 / 2 + l 



;2 + v^)e^-(l + v^) 



Returning now to the general case, we notice that the selfdual manifold 
always contains two special points for which the behaviour of the M coupled 
models can be related to that of a single Potts model. At the first such point, 



q 



m/2 



for m = 0, 1, 



M 



one has K\ = log(l + y/q) and K m = for m = 
models simply decouple. The other point 

b m = S(m, 0) for m — 0,1, 



2,3, 



(13) 



M, whence the M 



M 

T 



(14) 



corresponds to K m = for m = 1, 2, . . . , M — 1 and K M = log(l + q 
whence the resulting model is equivalent to a single g M -state Potts model. 
Evidently, for M = 1 these two special points coincide. 

Specialising now to the case of a regular two-dimensional lattice, it is well- 
known that at the two special points the model undergoes a phase transition, 
which is continuous if the effective number of states (q or q as the case may 
be) is < 4 fl~5|| . In Ref. || the question was raised whether one in general 
can identify further non-trivial critical theories on the selfdual manifolds. 
In particular it was argued that for M = 3 there is indeed such a point, 



7 



supposedly corresponding to a conformal field theory with an extended 5*3 
symmetry. 

To get an indication whether such results can be expected to generalise 
also to higher values of M, we have numerically computed the effective central 
charge of M = 4 coupled models along the two-dimensional selfdual surface. 
We were able to diagonalise the transfer matrix for strips of width L = 
4, 6, 8 lattice constants in the equivalent loop model. Technical details of the 
simulations have been reported in Ref. ||. Relating the specific free energy 
fo{L) to the leading eigenvalue of the transfer matrix in the standard way, 
two estimates of the effective central charge, c(4, 6) and c(6,8), were then 
obtained by fitting data for two consecutive strip widths according to 

/o(£ Wo(oo) - ^ + • • • . (15) 

A contour plot of c(6, 8), based on a grid of 21 x 21 parameter values for 
(bi, 62), is shown in Fig. [l]. The data for c(4, 6) look qualitatively similar, but 
are less accurate due to finite-size effects. We should stress that even though 
the absolute values of c(6, 8) are some 4 % below what one would expect in 
the L — > 00 limit, the variations in c are supposed to be reproduced much 
more accurately [§]. On the figure q = 3, but other values of q in the range 
2 < q < 4 lead to similar results. 

According to Zamolodchikov's c-theorem fL0| , a system initially in the 
vicinity of the decoupled fixed point (61,62) = (^,<l), shown as an asterisk 
on the figure, will start flowing downhill in this central charge landscape. 
Fig. Q very clearly indicates that the decoupled fixed point acts as a saddle 
point, and there are thus only two possibilities for the direction of the initial 
flow. 

The first of these will take the system to the stable fixed point at the 
origin which trivially corresponds to one selfdual g 4 -state Potts model. For 
q = 3 this leads to the generation of a finite correlation length, consistent 
with c c fj = in the limit of an infinitely large system. As expected, the flow 
starts out in the 62 direction, meaning that it is the energy-energy coupling 
between layers (K2) rather than the spin-spin coupling within each layer (Jfi) 
that controls the initial flow. 

More interestingly, if the system is started out in the opposite dirrection 
(i.e., with K 2 slightly positive) it will flow towards a third non-trivial fixed 
point, for which the edge weights tend to infinity in some definite ratios. 
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Figure 1: Contour plot of the effective central charge c(6, 8) along the self- 
dual surface (bi, b 2 ) for four coupled three-state Potts models. The decoupled 
fixed point is shown as an asterisk, and renormalisation group flow lines are 
sketched as a guide to the eye. 
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[Exactly what these ratios are is difficult to estimate, given that the asymp- 
totic flow direction exhibits finite-size effects.] Seemingly, at this point the 
central charge is only slightly lower than at the decoupled fixed point, as 
predicted by the perturbative renormalisation group ||. From the numerical 
data we would estimate the drop in the central charge as roughly Ac = 0.01 
- 0.02, in good agreement with the perturbative treatment which predicts 
Ac = 0.0168 + 0(e 5 ) §. 

All of these facts are in agreement with the conjectures put forward in 
Ref. H, and in particular one would think that this third fixed point corre- 
sponds to a conformal field theory with a non-abelian extended S4 symmetry. 

Finally, the numerics for q = 2 (four coupled Ising models) is less con- 
clusive, and we cannot rule out the possibility of a more involved fixed point 
structure. In particular, a c = 2 theory is not only obtainable by decou- 
pling the four models, but also by a pairwise coupling into two mutually 
decoupled four-state Potts (or Ashkin- Teller) models. Indeed, a similar phe- 
nomenon has already been observed for the case of three coupled Ising models 

§■ 
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